Exact results for the Wigner transform phase space densities of a two-dimensional 
harmonically confined charged quantum gas subjected to a magnetic field 
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Closed form analytical expressions are obtained for the Wigner transform of the Bloch density matrix and for 
the Wigner phase space density of a two dimensional harmonically trapped charged quantum gas in a uniform 
magnetic field of arbitrary strength, at zero and nonzero temperatures. An exact analytic expression is also 
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occupied, is also examined, and useful approximate expressions for the spatial and momentum densities are 
given. 
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I. INTRODUCTION 

Considerable interest has been shown in the study of the 
properties of the so-called low dimensional systems. The ad- 
vances in nanotechnology allows nowdays the realisation of 
quasi-two dimensional systems like quantum dots iH In 
a different context, the experimental achievement of trapped 
ultra-cold atom gases allows to study quantum mechanical ef- 
fects of quantum statistics in such gases Isf]. The above men- 
tionned physical systems have originated a great volum of the- 
oretical work in order to understand such fascinating world in 
a reduced physical space |4]. In this context, using the canoni- 
cal Bloch density matrix as a tool, exact analytical expressions 
have been obtained for the particle and the kinetic energy den- 
sities in spatial coordinates at zero and nonzero temperatures 
lUt]. Very recently, this method has been generalized to take 
into account the effect of a uniform perpendicular magnetic 
field on a confined charged two dimensional quantum gas |@] . 
In the present work, we are interested in obtaining exact an- 
alytical expressions for the Wigner transforms of both, the 
canonical Bloch density matrix and the first-order density ma- 
trix. Our interest in the Wigner tranform is based on the fact 
that it provides a useful tool to study various properties of 
many -body systems |7]. Besides the fact that it allows a refor- 
mulation of quantum mechanics in terms of classical concepts, 
and it is also used to generate semi classical approximations 
i^l^, the Wigner transform may help to gain a better under- 
standing of the properties of the system. It is also very in- 
teresting because with the recent progress of the experimental 
techniques, phase space densities can be nowdays measured 
for certain quantum systems dF 



The canonical Bloch density matrix is defined as 
C(r,r', j3) = Jlj(Pjir)(pj{r')exp{-j5ej), where (pj{r) and ej 
are eigenfunctions and eigenvalues of a one particle Hamilto- 
nian H associated to the system. Here, )3 is to be interpreted as 
a mathematical variable, which in general, is taken to be com- 
plex, and not necesarily the inverse temperature. The Bloch 
density matrix is of particular interest since its knowledge 
enables the first-order density matrix p(r,r') to be found, 
through the inverse Laplace transform llllll . In fact at T = 0, 
the first-order density matrix, for a given Fermi energy X, is 
given by 



(1) 



The system we are going to study is a harmonically confined 
charged atom gas in a two dimensional x-y plane subjected to 
a perpendicular homogeneous magnetic field B ^ Bk, taken 
along the z axis. The one particle Hamiltonian is then given 
by 



H 



1 

2m* 



-V + -A 

i c 



(2) 



with r = X +.y , A = j(B x r) is the vector potential, m* 
and — e are respectively the effective mass and the charge of 
the particle and a)o is the oscillation frequency of the confining 
potential. For the Hamiltonian under study, a closed analytical 
expression was obtained long time ago for the corresponding 
Bloch density lfl2ll . Here we rewrite it in the following useful 
form 



C(r,r',)3): 



m*a/2nh 
sinhj3/in 



exp 



-2m* a/ h 

sinhjS/ifl 



E} sinh — sinh 



-cosh 



cosh ■ 



— sinhjSricOi 



(3) 



2 



where R= (r + r')/2 and s = r — r' are, respectively, the 
center of mass and relative coordinates, and 

(^L = ^M= \fi[+^^,a± = a±coL. (4) 

COl is the Larmor frequency and Q.± are two frequencies 
that correspond to excitations in the center of mass motion- 
the so-called "Kohn modes" ifTsll . Note that, the Hamil- 
tonian in Eq. (2) has the same partition function, Z = 
1/ [4sinh()3/!n+/2)sinh(j3/!f2_/2)], as that an anisotropic 
two dimensional harmonic oscillator with frequencies Q. and 
Q.+. 

The rest of the paper is organized as follows. In the next 
section, we calculate the Wigner transform of the Bloch den- 
sity, and alternative useful analytical forms for such Wigner 
transfom are also derived. The Wigner phase space density 
matrix is calculated at zero and nonzero temperatures in sec- 
tion 3, showing some interesting plots of it. In section 4, we 
derive a closed analytical form, in terms of Laguerre poly- 
nomials, for the so called autocorrelation function. The high 
magnetic field strength case is examined in section 4. In the 
last section, a summary and outlook are given. 

11. THE WIGNER TRANSFORM OF THE BLOCH 
DENSITY MATRIX 

In the following we shall calculate the Wigner transform 
of the Bloch density matrix given in Eq. (3). The Wigner 



transform of an arbitrary one particle operator A, defined by 
its matrix elements in spatial coordinates A{r + |,r — |), is 
the following function Avi/ of the phase space variables r and 



Aw(r,p)= / A{r + s/2,r-s/2)e-'P-''^^' ds, (5) 



being its inverse transform 



A(r + s/2,r - s/2) = Aw{r,p)e'P-^/'. (6) 



According to Eq. (6), the local part of the operator A can be 
computed as 



Air,r)^Air)^l^A^ir,p). (7) 



We can now calculate, by making use of Eq. (5), the Wigner 
transform of the Bloch density matrix (3). Let us take the 
Wigner transform of C(r + |,r— |,j3) and callitCw(^',P,j3), 
so that 



m*a/2nh 
sinh^TiQ. 



exp 



2m*Q. 



exp 



hsiDhilihD.) 
2/jsinh(/3ria) 



sinh 



sinh - 



cosh 



2 



cosh 



2 



m* Hsinhp ficoi 
hsinhiphO.) 



(fcxr) + | 



The above two dimensional integral can be easily evaluated where for notational simpUcity we have introduced the fol- 
by using the well known identity lowing functions of j3 



ds e 



to obtain the result 



(9) 



2m* n sinh 



sinh 



TisinhiPhD.) 
sinh()3rii2) 



2m*hQ.cosh cosh ^^y— 
m*i2sinhj8/!«L 
smh{j3hQ.) 



(11) 



,,„, 1 ,„w , ,„,,, „7 It can be easily checked that when the magnetic field IS absent, 

p-f{P)r- „-g{P)ip+uil3){kxr)) , ^ r^ 

Cw{r,p,l5) — - - , (10) ^° "^^^ ®L = then i2+ = ii_ = tOo and i2 = coq, Eq. (10) 

cosh cosh yields to the correct Wigner transfom for a harmonic oscillator 



3 



in two dimensions, that is 



takes the following simple form 



exp 



2tanh 



h(0o 



p m CO, 

2m* 



V2 



cosh 



2 priiDQ 



(12) 

For the case of an unconfined system subjected to a magnetic 
field, i.e. coq ~ 0, then D. ~ col, ii- = 0, Q.+ = 2(0l, and Eq. 
(10) reduces to 



C^-'{r,p,P) 



exp 



tanhjSTicoL / ^ 
2m*Ti(£)L V c 



cosh pliCOL 



(13) 



which is the correct expression of the Wigner transform II15I1 . 



G(r,p,B) = tanh tanh . 

(16) 

Substitution of this result into Eq. (15), leads to the factorized 
analytical form 



Cwir,p,p) 



exp 



tanh 



ha 



cosh 



exp 



tanh 



n^iHo-^L, 



cosh 



(17) 



A. Alternative forms for the Wigner transform of the Bloch 
density matrix 

In the following, we present two alternative analytical 
forms of the result in Eq. (10), which can be rewritten as 

Cw(r-,p,/3)= ^^^ ^ exp[-G(r,p,/3)] (14) 

with 



cosh ^"2^ cosh 



G{r,p,l5)^{f + gu^)r^+gp^ + 2guL, 



(15) 



and L, is the component of the orbital angular momentum 
along the z axis. We show in Appendix A that G{r,p,j5) 



J 



where 



Ha 



P 

2m* 



m*Q? 



(18) 



Let us now obtain a third closed expression for Wigner trans- 
form of the Bloch density. For that purpose, we use the fol- 
lowing expansion in terms of Laguerre polynomials iflill 



exp(-xtanh,) ^^^_,£ ^^ ^^ 
coshy „fo 



-2yin + -] 



(19) 

for.jc = (HoztClL-) /TiQ. andy = l5TiQ.±/2, Eq. (17) becomes 



Cw(r-, p,j3) =4exp 



2Ho 

hD. 



^ ^ (-1)"+'"L„ 

n=Om=0 



f2{Ho + DL^\j ( 2{HQ-aL„ , 



(20) 



where 



£„,,„ = hD+{n + 1/2) + M2_ (m + 1 /2) 



(21) 



are the eingenvalues of the Hamiltonian (2). To our knowl- 
edge, the results given in Eqs. (10), (17) and (20) are new and 
seem not to have been reported before in the literature. 



m. QUANTUM WIGNER PHASE SPACE DISTRIBUTION 
AT ZERO AND NONZERO TEMPERATURES 



A. Phase space distribution at zero temperature 

Having established in the previous section various analyti- 
cal forms for the Wigner transform of the Bloch density, we 
shall now calculate analytically the expression for the quan- 
tum Wigner phase space distribution or Wigner transform den- 
sity of the first-order density matrix p(r,r'). Let Pw{f^p) 



denote such density, defined as 



pw(r,p)= p{r + -,r--)e-'P-^/'^ds, (22) 
Jr2 2 2 

The above distribution can also be obtained through the use of 
the Wigner phase space version of (1), that is 

2ni Jc-ico p 

Inserting Eq. (20) into Eq. (23), and performing the inverse 
Laplace transform IitIi 



,-r'°°^.'^(^-^"'")=0(A -£„,„,), (24) 
2ni jc-ioc p 

where is the Heaviside step function, we find 



4 



n=Oin=0 



2{Ho + Q.L, 
HQ. 



TiQ. 



(25) 



r 



Due to the presence of the step function, the quantum numbers 
n, m are restricted to TiQ.+ (n + 1 /2) + 7iQ_ {m+l/2) < X. The 
highest allowed value for n, A^+, is given by 



Int 



£27 



(26) 



where lnt{x) denotes the integer part of x > 0. For a given 
allowed value of n, the maximum allowed value of m, is 



where A is the Fermi energy. Taking the Wigner transform of 
Eq. (30), we get 



Pw(r,p)= £ ^[(/)„,,„(r + |)C„(r-|)] 0(1 -£„,,„). 

(31) 

Here the symbol W stands for Wigner transform. Comparing 
this result with Eq. (25), we deduce that 



= Int 



£2^ 
IT 



£2 
£27 



(27) 



Therefore, the density distribution in Eq. (25), can be rewrit- 
ten as 



<l>nAr + ^)<l>L,{r-^)]=4e-Tii-ir"' 



2«o ^+ ^- 



pw(r,p) = 4e-^^ ^ (-!)"+"%, 

(1=0 m=0 



2(//o + £2L-.) 

r7£2 



Thus, we have found the Wigner transform of the product 
0«,m('" + s/2)0,t„,(^~ •s/2) without the explicit use of the sin- 
gle particle wavefunctions. As stated before this result will 
immediately be used in the following subsection. 



Notice that the above density depends not only on the moduli 
|r I and but also on the relative angle 9 between r and p, 
i.e., pw{r,p) = pw{i',p, 9). For N ^20 particles, Fermi en- 
ergy X = 6.35TiC0l and COq/cOl = 1, we display in Fig. 1 (a) to 
(d) this Wigner phase space density for 9 = 0, ;r/6, ;r/3, 7r/2, 
respectively. 

In the absence of a magnetic field, the system becomes of 
a pure harmonic oscillator with frequency coq and thus all the 
£2 are all equal to coq, we show in Appendix B, that the above 
density has indeed the correct limit, given by lUSll 



B. Wigner phase space distribution at nonzero temperatures. 

Here, we shall generalize the result obtained in (25), valid 
for r = 0, to nonzero temperatures. We start with the defini- 
tion of the first-order density matrix p{r + s/2,i — s/2,T) at 
temperature T in terms of the normalized single particle wave 
functions 0„_„, , which reads for Fermions 



Pw='ir,p) 



-4e 



2Ho 



M 



p=0 



1 (^1^ 

" VlCOo 



(29) 



where is the generalized Laguerre polynomial of order one 
and the quantum number M is related to the Fermi energy by 
X=hcoo{M+l). 

In what follows, we shall deduce from Eq. (25) a result 
which will greatly simplify us the treatment when we will 
deal, in the next subsection, with the finite temperature case. 
Let (j)„j„{r) denote the eigenfunction of the Hamiltonian (2) 
with eigenvalues given in Eq. (21). In terms of the single 
particle wavefunctions, the first-order density matrix is 



(30) 



<^>„.m(r + f)(/)*,,„(r-f) 



n=Om=0 exp I ^YgT ^ ' + ^ 



where 



1 



-1 



(33) 



is the Fermi distribution function 



for the level energy £„.„,, kg is Boltzmann's constant and jj. the 
chemical potential. Taking the Wigner transform of both sides 
in Eq. (33) and using obvious notations, we get 



n=Om=0 



W[(^„,„{r + s/2)^*Jr^s/2)] 



exp 



1 



(34) 

where we have used the fact that, the Fermi distribution is not 
affected by the Wigner transformation. Substituting Eq. (32) 
into Eq. (34), one arrives at 



5 



£11.111 —{J 

e V 



(35) 



r 



Since, in the T ^ limit the Fermi distribution function be- 
comes the step function, with A = /i (r = 0) as the Fermi en- 
ergy, that is. 



1 



1 



■ ©(A —£■„,,„), 



(36) 



the result (35) reduces to the correct zero temperature limit 
given in Eq. (25). As can be seen in (35), the Fermi distribu- 
tion function enters in a simple way in the expression of the 
phase space distribution. This suggests to examine a similar 



situation for the case of bosons. In this case, the first-order 
density matrix in spatial coordinates at temperature T is 



0«,m(r- + f)0;,,„(r-f) 



n=Om=0 



(37) 

where we have included the Bose distribution function. Fol- 
lowing the same derivation as done for Fermions, one imme- 
diatly gets for the phase space density of bosons at finite tem- 
perature, the result 



J 



p^(r,p,r)=4.-^£ £(-l)"+'"L„ 



«=0m=0 



2{Hq + £IL, 

no. 



2(j/o -HLQ 
TiD. 



-1 



(38) 



This last equation may constitute a useful starting point to 
study thermodynamical properties in phase space of charged 
Bose gas, in particular at low temperatues. 



IV. THE AUTOCORRELATION FUNCTION 

The autocorrelation function, also called the reciprocal 
form factor lfl9ll . is known to provide information on the off- 
diagonal part of the density matrix p(r,r') and is defined in 
spatial coordinates as 



B{s) 



exp 



ip ■ s 



n{p)dp 



(39) 



with s ^ r — r' and n (p) is the density profile in momentum 
space. The latter is defined by a similar relation as in Eq. (30), 
where one has to convert the normalized spatial wavefunctions 
^„_,„ into their analogues ^„,,„(p) in momentum space, that is 

= E £ 0n,m(p)C,(p)e(A - E„.„,). (40) 

;i=0;n=0 

On the other hand the momentum density n{p) can also be 
obtained through the Wigner phase space distribution 



= / Pwir,p) 



dr 



(41) 



and is normalized to the total particle number of the system 



n{p)dp = N. (42) 

Therefore, it follows from Eq. (39), that B(0) = N. In the 
following we shall derive a closed analytical result for B{s). 
To do so, we first insert Eq. (23) into Eq. (41), to obtain 



"(^^-2k£r'^TL(^^'^^^'^'^^ '''' 



to carry out the r integration, we use the analytical form of 
Cii/(r,p, j3) given in Eq. (10) and we rewrite it as follows 



Cw(r,p,j3) 



cosh ^ cosh ^ 



exp 



X exp 



f + gu 
gu 



(44) 

2 



f + gu' 



(pxk] 



where we have used (p x k)^ = p^, since p is a planar vector. 
The above result can now be inserted into Eq. (43), to obtain 



6 



IniJc-ioo ^ (2;rri)2j3 



cosh£^cosh^(/ + ^«2: 



■ exp 



Inti 



2sinh 



sinh 



-exp 



where we have used (11). Putting this result with this present form into Eq. (39) and using Eq. (9), one then finds 



B(s^ = — r'°°dQ "^P^^^^ 

^ ' iTtiJc-io. 27thm*D.psir)h{l5hD.) 



dp exp — 



2 smh ~ + sinh - — 



1 

iTti 



dp 



exp(j3A) 



1 



\ hm*D.smh{lihD.) 

m*Q. sinh(j3fin) 



2^ 2 'P-S 



4 sinh 



sinh 



- exp 



„2 



%h sinh ^ Sinh ^ 



(45) 



(46) 



Remember that, Q.= (i2++i2_)/2, so that 



(50) by the Fermi function 



sinh(j3rii2) 



smh ^ T smh 



= coth 



coth 



(47) 



2 2 

plugging this result into the exponential of Eq. (46), to get 



B{s) 



1 

2m Jc 



dp- 



exp ( -^coth 



exp 



m'Q. 
8/i 



Coth^U2 



2 sinh 



2 sinh ^ 



(48) 



At this level, we can carry out explicitly the inverse Laplace 
transform by first using the following expansion in terms of 
Laguerre polynomials lfl6ll 



exp 



-a: coth 



sinh(^; 



2e-' £ L„(2x) e-i5f'"±("+i) (49) 

n=0 



and followed by Eq. (24), one then finds 



m,n=0 



m*Q. 



I L,r, 



m*Q. 



(50) 

Notice that, an interesting feature of the above autocorrelation 
function is that it is expressed in terms of Laguerre polynomi- 
als with same arguments, note also that it is isotropic in spatial 
coordinates, i.e, depends only on the length |s|of the vector s. 
Setting s = 0, one obtains B(0) = A^, as it is required. In Fig. 2 
we display the T ~Q spatial dependence of the autocorrelation 
function for N = 20 particles, choosing the following values 
of the parameters (Oq/cOl = 1 and Fermi energy X = 6.35h(OL. 

Finally, let us come to the finite temperature expression of 
B{s). In a similar way as was done in the previous section, 
one can immediatly write down its expression. For Fermions, 
this is simply achieved by replacing the &{X — E„^„,) in Eq. 





( m*n 

Ah 




K, 1 




exp 







(51) 



For Bosons, all that is required is the replacement of the Fermi 
function by the Bose function 



B{s)=e-^' ^ 



Ln 


( m'a 

Ah 








exp 







(52) 



V. STRONG MAGNETIC FIELD CASE 

The strong magnetic field case at low temperature is of 
particular interest in quantum dots. In this limit only a few 
Landau levels are occupied and the magnetic field length 
/ = (/ic/eB)'/^ is small leading to slowly varying confining 
external harmonic potential on the scale of /. In what fol- 
lows we shall examine the strong magnetic field (SB) case, 
where one has coo/ol <C 1 then H w ©l, i2+ w 2a)L and £1 w 
COq / (2(Bl) which yields respectively for the functions given in 
Eq. (11), to the leading order 



M(j3) Wm*(BL. 



(53) 

Substituting this results into Eq. (10), gives immediately the 
result for the Wigner transform of the Bloch density 



■exp 



tanhjS/icOL / e 
2m*h(x>L 



P 



(54) 

In this case, the net result we get is a product of the Bloch den- 
sity of free charged particles in magnetic field [see Eq. (13)] 
and an exponential factor limiting the spatial distribution. Let 
us now, calculate the corresponding Wigner phase space den- 
sity. Inserting Eq. (54) into Eq. (23) and making use of Eqs. 



7 



(19) and (24) , one obtains 



Pwir,P) = 2exp 



H, 



tnagn 
TiCOt. 



li^irul^J^) (55) 



H=0 



nTG), 



'0^2 



where = (p+ (e/c)A)^/(2m) is the Hamiltonian for a 

particle in the presence of the magnetic field alone. One im- 



mediately recognizes in the argument of the Heaviside func- 
tion the discrete Landau level energies (2n + I)/!©!- As 
can be seen, the phase space density in above has a sim- 
ple analytical form, therefore we can easily obtain, in this 
high magnetic field limit, the corresponding density matrix 
p(r + s/2,r — s/2)in spatial coordinates. To do so, we make 
use of the inverse Wigner transformation. According to Eq. 
(6), one has 



J 



p^^(r + s/2, r - s/2) = 2 52 (- 1 )"© ( A - (2« + 1 )ritOz, 



n=0 



(0, 



'Oj.2 



r 



dp 



tinuign 



2H, 



magn 



Ticol 



(56) 



The last integral can be carried out as follows. Denoting it by 
/ and using the canonical momentum K = p + {e/c)A, one 
obtains 



I — e 



K2 {27th) 



m*hcoi J ' 



(57) 



and changing to the variable t ~ K/^m*hcoi, we get 



(2;rri)2 



j^^J te '-TL„{t^)Jo{x/m*cOL/fits)dt, (58) 





where we have made use of the relation ll 



Jo 



to get the last line, 7o(j^) being the Bessel function. The fol- 
lowing relation 1II61I 

rco ^ 2 

xe-^Ln (x2) jQ{xy)dx = (-l)"e-^L„ (/) 
helps to perform the integral in Eq. (58), to find 

l=t}le-^lA.s^--^s^LJ"-^sA (59) 



27r/2 " " ~" V Ti 

Substituting this result into Eq. (56), yields 



2nf 



J7=0 



m COl 2 
s 



x0 (A - (2« + l)hcOL - m*(0^r^/2) . 
The local density is obtained by setting s = 0, 

P''(^) = 2^ L - (2« + - '^^') ■ (60) 



n=0 



In their study of a two dimensional electron gas subjected 
to a magnetic field and partially confined by a harmonic po- 
tential, the authors of Ref. 12011 obtained a similar result for 
the spatial density using a different approach (remark that 
their parabolic potential is taken only in the x direction, i.e, 
V{x,y) = m* (OqX^ /2). As noticed by these authors, the den- 
sity contains compressible and incompressible regions. 

In order to rewrite Eq. (60) in a more compact form, we use 
the following identity relating the Heaviside and the integer 
part functions 



Then, p(r) becomes 
1 



-«) = 0(x+l)lnt(x+l) 



(61) 



p(r) 



27r/2 







xint 



2Ti(0l 

A — m*(olr^/2 + Ticol \ 



2h(0L 



(62) 



In Fig. 3 we plot the above zero-temperature spatial density 
for the case of N = 200 particles with parameters coq/col = 
0.2048 and Fermi energy A = 4.1JiC0l 

For ultra strong magnetic field, such that all the particles 
reside in the lowest Landau level (LLL), Eq. (60) reduces to 

P^^^(r-) = ^0 (A - Ticol - m*&lr'l2) (63) 

Before closing this section, it is intersting to calculate the 
momentum density or the momentum distribution «(p) in 
the strong magnetic field case. This important distribution 
was already introduced in the previous section but its calcu- 
lation has not been fully completed since this density was 
used there as an intermediate to obtain the autocorrelation 
function. In a short, we start from its expression given in 
Eq. (46) and introduce the strong magnetic field approxi- 
mations we used above, namely sinh (j3/!i2+/2) w sinh jS/iOi, 
sinh(j8r!i2_/2) w j5h(0^ / {4col) and sinh()3rii2) w sinhjihcoL 



8 



to get 



dp 



exp(j3A) 



{27th)m* coip sinhjS/zOL 



(64) 

At this level, it is easy to perform the inverse Laplace trans- 
form by using the expansion 

■ . , = 2 £ exp [-(2« + 1)I5Ti(0l] 



phase space density at zero and nonzero temperature. Our re- 
sults are valid for arbitrary magnetic field strengths and hold 
for both Fermions and Bosons. For the system under study, we 
have found simple and exact analytical expression for the so 
called autocorelation function. The high magnetic field case 
has been examined. Our investigation in phase space comple- 
ment the recent works in spatial coordinates. The results we 
obtained would constitute useful starting point for the study, 
in phase space, of thermodynamical properties in the field of 
cold atom gases. 



with Eq. (24), which results in 



„SBf 



21^ 



n=0 



p) = — 2 ^ ( A - i2n + l)hC0L- 



'0 



2m* (0\ 



(65) 



Like the spatial density in Eq. (60), the momentum density ex- 
hibites the same structure consisting of a series of wide steps 
in momentum space. 

The above results have been obtained for spinless charged 
particles at T = and the generalization to finite temperatures 
can be done without any particular difficulties. 
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VI. SUMMARY AND OUTLOOK 

We have derived some simple exact closed expressions for 
the Wigner transform of the canonical Bloch density of two- 
dimensional harmonic oscillator in a uniform magnetic field. 
We have also obtained exact analytical form for the Wigner 



APPENDIX A 

The purpose of this appendix is to derive the expression of 
the function G(r, p,j3) given in Eq. (17). First, we evalu- 
ate the various functions of j3 in Eq. (15), namely f/()3) = 
(/(j3)+g(j3)«2(j3)), g(j3) and y(j3) = 2g(j3)M(j3). Using 
Eq. (11), we getfor;7(j3) 



J 



2,„*£2sinh^sinh^ 

U{B) = 2_ 2 

risinhCfl 



m*i2sinh^j3/!a)L 



Z Z I ^ ; 

&mh{pha) 2/jsinh(j3/jn) cosh ^ cosh ^ ^ sinh/3r7i2 



m*a/2fi sinh()3rii2_)sinh(j37!i2+) + sinh2j8fia)L 



. ana- , am^ 
cosh — cosh i-T-i 



r 



iAl) 



From sinh(j37!£2_)sinh(j3rii2+) = \ (cosh(j37!(£2_ - 



cosh()37i(i2_ and Q. +n 

—2(x)l, we obtain 



2a , a - n+ 

m*Q. [cosh2j3/!n - coshj3fi©,. + 2 sinh^ PHcOl] 



4risinh(j3rii2)cosh cosh^^ 
m*Q.smhj3nQ. 



2h cosh ^^'P cosh 



2 2 

If we use again Q. + D.+ = 2£2, the above result becomes 



tanh ■ 



(A2) 



We now turn to the two remaining functions. For convenience, 
we rewrite the function g(j8), given in Eq. (11), as 



eiB) = tanh^- 



tanh 



(A3) 



For the function V{j3), one simply gets 

sinh2j3/!a)L 



/icosh 



cosh 



2 2 

and using i2+ — D. = 2coi, one ends with 
V{I5) = - f tanh ^iyi - tanh 



(A4) 



Substituting Eqs. (A1)-(A3) into Eq. (16), simple manipula- 
tions yield to the desired result (17). 



APPENDIX B 

In this appendix we shall show that, in the absence of mag- 
netic field, the phase space density in Eq. (28) reduce to the 
result in Eq. (29). In this limit Eq. (28) becomes 
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Here, = Int( ^ ) — 1 and A'- = Int ^ ^ — « — 1 ^ = — n. The physical meaning of is only but the quantum number of 

the last occupied harmonic oscillator shell, denoted by M in Eq. (29). The Hamiltonian Hq appearing in Eq. (Bl) refers to Eq. 
(18) but with Q. = C0(). Putting p = n + m, Eq. (Bl) rewrites 



and using the identity Y!m=o ^p-m {^) {y) (-^ + y) (see lll6ll ). the result in Eq. (29) is then recovered 
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FIG. 1 : Figures (a) to (d) correspond to plots of the Wigner phase 
space density pjy(r,p) = pjy(r, 6) for N = 20 particles at = 
0, ;r/6, ;r/3, ;r/2, respectively. We have chosen parameters coq/coi = 
1 with Fermi energy A = 6.35/iCOl. Lengths are plotted in units of 
the magnetic length I = yJhcjeB for r and in units of for the 
momentum p. 
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FIG. 2: The r = autocorrelation function B{s) = B{s) for A' = 20 
particles, with o\)/col = 1 and Fermi energy A = 6.35/;©^. Lengths 
are plotted in units of the magnetic length / = sJhcjeB. 
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FIG. 3: Plot of the spatial density p(r) = p(r) in units of ilnfy'^ 
at r = in a harmonic oscillator potential with coq/coi = 0.2048 
and Fermi energy A = AAhcoi, corresponding to A' = 200 particles. 
Lengths are plotted in units of the magnetic length / = sJhcjeB. 



